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Abstract. We theoretically investigate various beyond mean- field effects on Bose 
gases at zero temperature featuring the anisotropic and long-range dipole-dipole 
interaction in addition to the isotropic and short-range contact interaction. Within 
the realm of the Bogoliubov-de Gennes theory, we consider static properties and low- 
lying excitations of both homogeneous and harmonically trapped dipolar bosonic gases. 
For the homogeneous system, the condensate depletion, the ground-state energy, the 
equation of state, and the speed of sound are discussed in detail. Making use of the local 
density approximation, we extend these results in order to study the properties of a 
dipolar Bose gas in a harmonic trap and in the regime of large particle numbers. After 
deriving the equations of motion for the general case of a triaxial trap, we analyze 
the influence of quantum fluctuations on important properties of the gas, such as 
the equilibrium configuration and the low-lying excitations in the case of a cylinder- 
symmetric trap. In addition to the monopole and quadrupole oscillation modes, we 
also discuss the radial quadrupole mode. We find that the latter acquires a quantum 
correction exclusively due to the dipole-dipole interaction. As a result, we identify 
the radial quadrupole as a reasonably accessible source for the signature of dipolar 
many-body effects and stress the enhancing character that dipolar interactions have 
for quantum fluctuations in the other oscillation modes. 
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1. Introduction 

The experimental realization of Bose-Einstein condensation in a ^^Cr sample by the 
group of Tilman Pfau in 2005 triggered much experimental and theoretical work in the 
field of dipolar quantum gases [1]. Chromium atoms possess magnetic moments of 6 
Bohr magnetons (/ib) so that the anisotropic and long-range dipole-dipole interaction 
(DDI) between them is 36 times larger than the one between alkali atoms. Therefore, 
taking the influence of the DDI into account besides that of the isotropic and short-range 
contact interaction is essential for the correct physical description of chromium Bose- 
Einstein condensates (BECs). Up to date a few experimental signatures of the DDI in 
BECs have been identified. The most striking ones have been found in chromium such 
as the modified time-of-fiight dynamics [2], the strongly dipolar nature of a quantum 
ferrofiuid [3], the d-wave Bose nova explosion [1], and the modified low-lying excitations 
[5]. In addition, the DDI has also been observed in rubidium [6] and in lithium [7] 
samples. For recent reviews on the physics of dipolar Bose gases see Refs. P [HI ITO] . 

In the case of the chromium experiments, careful comparisons with mean-field 
theory [HI [12] were carried out and brought remarkable agreement as a result. It 
is important to remark that the dipolar mean-field theory is based on the construction 
of the corresponding pseudo-potential [131 [H]- 

Among the strongly magnetic atoms, an important place is occupied by dysprosium, 
which possesses the unsurpassed magnetic dipole moment of 10 /ie- Recently, two major 
experimental achievements could be obtained with the trapping [15] and the subsequent 
Bose-Einstein condensation [16] of ^^^Dy by the group of Benjamin Lev. At present, 
Feshbach resonances are being searched for in dysprosium, which would provide a tuning 
knob for the relative strength of the DDI with respect to the contact interaction. It is 
also worth mentioning the possibility of using erbium as a strong magnetic atom. With a 
magnetic dipole moment of T/^b and a mass of 164 atomic mass units, erbium represents 
a promising candidate for studying dipolar physics [1^ . 

However, atomic systems are not the only dipolar quantum systems under current 
investigation. Indeed, the recent successes in producing and cooling heteronuclear polar 
molecules down to their rovibrational ground state by means of stimulated raman 
adiabatic passage (STIRAP) [TD], specially the cooling of KRb molecules [H] and 
the manipulation of their internal degrees of freedom [19] let us hope that quantum 
degenerate heteronuclear molecular systems will soon be available experimentally. And 
this is not all. By means of applied electric fields, lab-frame electric dipole moments 
can be induced in these molecules, thereby tuning the electric DDI over various orders 
of magnitude [20]. Typically, the DDI in polar molecules can be up to 10^ times larger 
than in atomic systems. 

As a natural consequence of the many important experimental successes, much 
theoretical effort has been dedicated recently to the investigation of strongly dipolar 
quantum gases. In the case of Fermions, one should mention at least the studies involving 
zero sound [211 122] and the Berezinskii-Kosterlitz-Thouless phase transition [23] in 
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homogeneous systems, superfluidity in trapped gases [21], and Wigner crystallization 
in rotating two-dimensional ones [25]. Particularly interesting are the investigations of 
fractional quantum Hall states [26l |27] . There have also been recent important studies 
involving bosonic dipoles considering, for example, finite-temperature effects [28], exotic 
density profiles [29], and the possibility of spin-orbit coupling [30]. Moreover, it was 
found that loading the system into an optical lattice leads to novel quantum phases for 
both bosonic [311 132] and Fermionic [33] dipoles. In the meantime, the first experiments 
with chromium loaded into an optical lattice have been carried out [M] . 

In view of the wide-range tunability of the DDI, polar molecules offer the possibility 
of testing dipolar mean-field theories all over and beyond their range of validity. For 
this reason, it is important to analyze theoretically dipolar systems beyond the mean- 
field approximation. Recently, we have briefiy discussed the infiuence of quantum 
fluctuations on trapped dipolar Bose gases [35]. In the present publication, we present 
in detafl the corresponding Bogoliubov-de Gennes (BdG) theory and apply it to both 
homogeneous and harmonically trapped gases, thereby emphasizing the importance of 
quantum fluctuations in strongly dipolar systems. 

This paper is organized as follows. In Section O we shortly discuss the BdG theory 
of a Bose gas at zero temperature containing a large number of polarized point 
dipoles. Section [3] is dedicated to homogeneous dipolar Bose gases, where we solve, at 
flrst, the Bogoliubov equations algebraically. Then we use this solution to study key 
properties of the system such as the condensate depletion and the beyond mean-fleld 
speed of sound. In Section IH we concentrate ourselves on harmonically trapped systems. 
By means of the local density approximation (LDA), we solve the BdG equations and 
derive the dependence of the condensate depletion and the ground-state energy on the 
system quantities such as the relative dipolar interaction strength and the contact gas 
parameter. With these results at hand, we work out a variational approach to superfluid 
dipolar hydrodynamics, which allows for deriving equations of motion for the Thomas- 
Fermi radii of the gas in the case of a triaxial harmonic trap. Section [5] is specialized 
to the case of a cylinder-symmetric trap and contains the solution of the equations 
of motion as well as discussions about beyond mean-fleld effects on both the static 
properties and the hydrodynamic excitations. Finally, in Section [6l we summarize our 
flndings and present the conclusions and perspectives of this work in view of future 
experiments. 

2. Bogoliubov-de Gennes Theory For Large Particle Numbers 

In this section we briefly present the most striking aspects of the BdG theory, which 
shall be applied to dipolar Bose gases in the following sections. Thereby, we emphasize 
the peculiarities which come about due to the non-local and anisotropic character of the 
DDI. 
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2.1. General Formalism 

Consider a gas of bosonic particles with mass M possessing a finite dipole moment 
at zero temperature. For definiteness, we consider the dipoles to be ahgned along the 
z axis. In this case, the interaction potential has a contact component V5(x) = 5'5(x), 
with the coupling constant g being related to the s-wave scattering length as through 
g = Airfi^as/M, and a DDI component which reads 

In the case of magnetic dipoles, the dipolar interaction strength Cdd is characterized by 
Cdd = /^o'^^) with fiQ being the magnetic permeability in vacuum and m the magnetic 
dipole moment, whereas for electric dipoles we have Caa = 4:7Td^ with the electric dipole 
moment d expressed in Debyes. As a whole, it is convenient to write the resulting 
interaction potential as 

3edd A ^ -2^ 



Mnt(x) = g 



(2) 



47r|xP \ |xP 

with edd = Cdd/'ig denoting the relative interaction strength. It is also convenient to 
introduce the dipolar length a^d = CddM/12'n'h'^ as a measure of the absolute dipolar 
strength, so that the relative interaction strength reads edd = fldd/os- 

To study the dipolar system within the BdG theory, we consider the total 
Hamiltonian H = Hq + -f^mt, which consists of a free and an interaction contribution. 
In general, the non-interacting part contains the kinetic and the trapping energy 

^0= frx^^(x)/io(x)^(x), (3) 



where \^"''(x) and ^(x) denote the usual bosonic creation and annihilation operators, 
respectively, and we have introduced the abbreviation 

ho{^) = -^ + UM. (4) 
Moreover, the interaction is included through 

Hint = \j^'^ jd'x'^\^)^\^')V,^t (x - x') ^(x')^(x), (5) 

where Mnt(x) is given explicitly by Eq. (jj]). We implement the BdG theory by 
diagonalizing the grand-canonical Hamiltonian H' = H — fiN, with the number operator 
= Jd^x\l/^(x)\l/(x) and the chemical potential /i. This is done by means of the 
Bogoliubov prescription \l/(x) = \l/(x) -\- 6ip (x.) , where the classical field ^'(x) represents 
the number A^o of condensate particles via A'^g = Jd^a;^'^(x)^(x), and the operator 
6tlj{x) accounts for the quantum fluctuations. 

By inserting the Bogoliubov prescription into the grand-canonical Hamiltonian, one 
can separate the contribution of the condensate and the one of the quantum fluctuations 
order by order in the fluctuation operator 6ip{x). Restricting the expansion to the zeroth 
order leads to the Gross-Pitaevskii equation 

vl/(x)/i= [/io(x) + (7|vl/(x)|2 + $dd(x)]M/(x), (6) 
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where the dipolar mean-field potential reads 

$dd(x)= /"dVVdd(x-xO|vl>(xOp. (7) 



The Gross-Pitaevskii equation (jH]) is the main tool in order to investigate mean-field 
properties of BECs. Although it can only be applied at very low temperatures and weak 
interactions, it has up to date been able to account for all experimental results obtained 
with dipolar BECs. 

As the interaction becomes stronger, one has to include the effects of quantum 
fluctuations. In order to do so, one flrst carries out the expansion of the grand-canonical 
Hamiltonian up to the second order in the fluctuations. Then, following de Gennes |36j . 
one introduces the expansion of the quantum fluctuations 

5^(x) = [K{^)a, + V:(x)dt] , (8) 

where the creation and annihilation operators dj, and also satisfy bosonic commuta- 
tion relations and the Bogoliubov modes are denoted through the index v. Here, it is 
important to exclude the ground state |0), which is deflned by djy|0) = 0, from the sum. 
This is denoted by the prime after the summation sign. Moreover, it is worth remarking 
that this expansion represents a canonical transformation if the Bogoliubov amplitudes 
Wj,(x) and Vj,(x) satisfy the condition 

jd'x [Z^:(xR,(x) - V:(x)V,,(x)] = 5,y, (9) 

which we shall, therefore, impose. Then, the resulting Hamiltonian will be diagonal, if 
the functions satisfy the BdG equations 

[e,-ifpi(x)]W,(x)= y"d32/Mnt(x-y) [M/(y)^(x)V,(y) + ^*(y)M/(xR(y)], 
-[£, + ifFi(x)] V,(x) = /d^yMnt (x-y) [vl/*(y)M/*(xR(y) + M/(y)M/*(x)V,(y)] , (10) 



where denotes the excitation energy of mode v. In the equations above, we have 
introduced the deflnition of the fluctuation Hamiltonian density according to 



i7Fi(x) = /io(x) 



/i + y"d3?/**(y)yi,t(x-y)vl/(y). (11) 



After diagonalizing the Hamiltonian, one can determine the number of particles in 
the many-body ground state |0) via N = (0|iV|0), which decomposes according to 

N = No + Y! /d=^a;V:(x)V,(x). (12) 

Thus, the total number of particles is a sum of the condensed and excited particles. The 
latter are moved from the one-particle ground-state to one-particle excited states due 
to the interaction, thereby depleting the condensate. 
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It should be also noted, that the ground state obtained under consideration of the 
quantum fluctuations differs from that of the Gross-Pitaevskii theory. More precisely, 
the ground-state energy is shifted by 

= - Jd'x [W:(x)i7Fi(xR(x) - V:(x)i/Fi(x)V,(x)] 

- |dVy"d3xl^,t(x-x')vl/*(x')vl/(x)[W:(xR(x')-V:(x)V,(x')]}. (13) 

Later one, we will use this correction to the ground-state energy f[T5]) as the starting 
point to determine the effects of quantum fluctuations upon BECs. 

In the general form presented here, the BdG theory is difficult to apply as the BdG 
equations are complicated to solve even numerically for a dipolar Bose gas due to the 
non- locality of the DDI [37]. Nevertheless, it is possible to find analytic approximate 
solutions in cases of special experimental interest. 

2.2. Thomas-Fermi Regime 

The BdG equations can be used to investigate the excitations of a Bose gas all the way 
from the harmonic oscillator regime, where interactions play no role, up to the Thomas- 
Fermi regime, where the interaction energy is much larger than the kinetic energy. In 
this paper, we are interested in the latter regime of a large number of particles and strong 
interactions, where the kinetic energy of the condensate can be neglected in comparison 
with the interaction and trapping energies. Thus, inside the condensed region, the 
time-independent Gross-Pitaevskii equation ([6]) assumes the following form 

/i = f/t,(x) + gnoi^) + $dd(x), (14) 

where the condensate density is identifled according to no(x) = \E'*(x)\l/(x). Indeed, 
the BdG Eqs. ffTOl) must be considered separately inside and outside the condensate. 
However, the solution for the external region implies V^(x) = 0. As a consequence of 
that, both the depletion and the correction to the ground-state energy vanish identically 
in this region. For this reason, we restrict our study to the condensate region. Under 
these circumstances, the BdG equations (ITUl) reduce with (IT^ to 

W,(x) = Jd'yV,^, (x-y) [vI/(y)vl/(x)V,(y) + vI/*(y)v^(x)W,(y)] , 

V,(x) = y"d=^|/l^„,(x-y)[v^*(y)vP*(x)W,(y) + vl>(y)vl>*(x)V,(y)]. (15) 

In the following sections, we will solve this set of coupled equations analytically 
for the case of a homogeneous dipolar Bose gas and that of harmonically trapped gas 
within the semiclassical approximation. 

3. Homogeneous Dipolar Bose Gases 

Even though homogeneous cold atomic systems cannot be realized experimentally, their 
study is of large importance. The reason for this is that these systems serve as a 



2 



2M 



2M 
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prototype for the experimentally relevant trapped cold atomic systems and often lead 
to the correct physical intuition with respect to their properties. Therefore, we start 
the application of the BdG theory by considering that the gas is enclosed in a volume 
V and that the field ^ is independent of position x. Therefore, the mean-field value of 
the chemical potential is given according to Eq. by 

/i = no lim V^nt (k) , (16) 

where the Fourier transformed of the interaction potential ([21) is written as [38] 

V^nt(k) = (7[l + edd(3cos2^-l)] (17) 

with 9 being the angle between the vector k and the polarization direction. Thus, the 
chemical potential is not uniquely defined due to the anisotropy of the DDI in three 
spatial dimensions. We will discuss how to solve this ambiguity in the following. For 
the moment, let us remark that, for a quasi two-dimensional dipolar system with the 
dipoles oriented perpendicular to the plane, the dipolar potential has a well defined 
value at the origin k = [39] . 

3.1. Bogoliubov Spectrum and Amplitudes 

Due to translation invariance, momentum is a good quantum number, so the excitations 
can be labeled with the wavevector k. In this case, the Bogoliubov equations (fT5|) turn 
out to be algebraic in Fourier space and read 

/l^k^ 

SkKu = ^TTT^^ + ^oMnt (k) [Wk + Vk] , 
ZM 

-£kVk = I^Vk + noVi^t (k) [Wk + Vk] . (18) 
Suitable algebraic manipulations allow to solve for both the Bogoliubov amplitudes 



2ek 



2M 

and the Bogoliubov spectrum [lO 



?^oMnt (k) - Ek 



(19) 



(20) 



Notice that, due to the condition (Q, the amplitudes are completely characterized by 
Eq. 

The Bogoliubov spectrum allows immediately to study the low-momenta properties 
of the system. Indeed, the sound velocity can be obtained by taking the limit k — t- 
of the spectrum. Again, the anisotropy of the DDI renders the limit dependent on 
the direction of the vector k, as its modulus tends to zero. For this reason, the sound 
velocity acquires a dependence on the propagation direction, which is fixed by the angle 
6 between the propagation direction and the dipolar orientation, and reads 



c{e) = jyVl + edd (3cos2^-l). (21) 
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Figure 1. Functions 23(2;) (red) and Q^ix) (blue), which govern the dependence of 
the condensate depletion and the ground-state energy correction on the relative dipolar 
interaction strength x = edd, respectively. 

This result represents the physics of the k = 0-mode as obtained from the Bogohubov 
theory. At this point, it is important to remark that, for edd > 1, the sound velocity 
may become imaginary depending on the direction of propagation. This instability of 
the system is an important characteristic of dipolar Bose gases which resembles the case 
of isotropic systems with attractive interactions |10]. 




3.2. Condensate Depletion 

Let us now study the number of particles expelled from the ground state by the inter- 
actions, i.e., the condensate depletion. As we are concerned with the thermodynamic 
limit, the quantum numbers become continuous variables and the summations can be 
replaced by integrals according to the prescription jH] 

Under these conditions, one finds that the condensate depletion is proportional to the 
square root of the gas parameter raa^ and reads [12] 

N - No 8 . 3x1/2 ^ . . , . 

= ^'''''^ ^ ^ 

The contribution of the DDI is expressed by the function Q3{x), which, for < a; < 1, 
is the special case / = 3 of 

a.(.) = (l-.)'.F, (24) 

Here, 2-^1(^5 7; ^) represents the hypergeometric function |13]. In Fig. [1] we plot 
the functions Qi{x) for I = 3 and I = 5 against x. It is worth remarking that the 
functions Qsi^x) and Q5{x) become imaginary for x > 1. Indeed, for edd > 1 the 
dipolar interaction, which is partially attractive, dominates over the repulsive contact 
interaction leading to the collapse of the condensate [1]. 
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For a condensate with pure contact interaction, this quantum depletion has never 
been observed due to difficulties in measuring the condensate density with sufficient 
accuracy and at low enough temperature. Including the dipole-dipole interaction only 
slightly increases the condensate depletion as given by formula (123|1 . Indeed, for the 
maximal relative interaction strength edd ~ 1 the condensate depletion can be about 
30% larger than in the case of pure contact interaction, so that the most important 
quantity remains the s-wave scattering length a^. Nonetheless, the establishment of 
Eq. (1231) is an important result from the theoretical point of view as it clarifies how the 
condensate depletion depends on the relative interaction strength edd- 

3. 3. Ground-State Energy and Equation of State 

Correspondingly, the presence of quantum fluctuations also leads to a correction of the 
ground-state energy of a dipolar Bose gas. Calculating this directly from the Bogoliubov 
theory yields initially an ultraviolet divergent integral. However, this can be repaired 
by calculating the scattering amplitude at low momenta up to second order in the 
scattering potential \4nt (k) and performing a proper renormalization procedure |44j . 
The final result can be obtained with the help of Eq. f|T3|l and reads 



with the auxiliary function Q5(edd) describing the dipolar enhancement of the correction, 
see Fig. [TJ Notice from Fig. [T]that Q^{x) varies from Q5(0) = 1 up to 25(1) ~ 2.60, 
so that the effect of the DDI is more significant for the energy correction f l25|) than 
for the condensate depletion fl2^ and offers, therefore, better chances for experimental 
observation. 

By differentiating the energy correction (I25p with respect to the particle number, 
one obtains the beyond mean-field equation of state 



In the case of a Bose gas with pure contact interaction, i.e., for edd = 0, this equation 
reduces to the seminal Lee- Huang- Yang quantum corrected equation of state |15]. 

It is worth remarking that, while the leading term in Eq. f l2B]l is anisotropic, in the 
sense that its value depends on the direction in which the limit k — )• is carried out, the 
sub-leading contribution from the quantum fluctuations is isotropic. The reason for this 
is as follows. As it accounts for the condensate, the leading term is evaluated at k — >■ 
and is, therefore, subject to the anisotropy in V^int (k) which is peculiar to the DDI. The 
second term, however, accounts for excitations with all k 7^ wavevectors. Therefore, 
it contains an integral over all the k-modes, which removes any possible dependence on 
the momentum direction. 




(25) 




(26) 
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Figure 2. Comparison between the sound velocities in the mean- field approximation 
([2T|) . plotted in dashed curves, and its quantum corrected version (|29)) . represented 
here in solid curves. The red and the blue curves are for e^d = 0.6 and edd = 1, 
respectively. 



3.4- Dipolar Superfluid Hydrodynamics 

Now that we have calculated the quantum correction to the equation of state, the 
corresponding correction to the sound velocity can be obtained by linearizing the 
superfluid hydrodynamic equations |16] around the equilibrium configuration. Consider, 
to this end, the continuity equation 

^^^ + V-Kx,t)v(x,t)] = 0, (27) 
where v(x, t) is the velocity field, together with the Euler equation 



dt 



yv(x,t)2 + /i(n(x,t)) 



(28) 



Linearizing these equations according to ?7,(x, t) = + 6n{:K,t) and assuming that the 
density oscillations have the plane wave form Sn{x,t) oc e«(kx-nt)^ ^^^^ obtains the 
corrected sound velocity 



fa\ /I , fi 2Z) n , 16v/a3nQ5(6dd) 
c{9} = C5\l I + edd [3 cos^ 9 - 1) -\ ^= , (29) 



with cs = \/ gn/M . Notice that (l29l) represents the extension of the Beliaev result for 
the sound velocity of Bose gases with short-range interaction [17] in order to include 
the DDI. The Beliaev sound velocity is usually displayed by expanding the square root 
with edd = in powers of the gas parameter a^n. Here, we prefer the form fl29l) because 
it is well defined for all directions and for all values of the relative interaction strength 
satisfying edd < 1- 

Notice that, in order to recover the mean-field sound- velocity (12T!) as it was given 
within the Bogoliubov theory, the limit in Eq. (l26l) was taken along the direction of 
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the sound propagation. Indeed, the hnearized hydrodynamic equations do capture the 
low-momenta physics and must, therefore, match the Bogohubov result at k — ?■ 0. 

Let us discuss the sound propagation for typical experimental values of the gas 
parameter a^n 10""^ of dipolar systems such as chromium [3]. The sound velocity as 
a function of the angle between the propagating wave and the dipole axis is plotted in 
units of cs in Fig. [2] for edd = 0.6 (red) and edd = 1 (blue). The continuous and dashed 
curves denote the quantum corrected and the mean-field velocities, respectively. The 
mean-field sound velocity for edd = 1 vanishes at ^ = tx/2. This is the signature of the 
instability of the system, as the partially attractive dipolar interaction dominates over 
the repulsive contact interaction for edd > 1- Including quantum fluctuations renders the 
sound velocity non- vanishing at this value of the interaction strength and propagation 
angle. The stability limit remains, however, unaltered due to the fact that the function 
Q^i^Ad) becomes imaginary for edd > 1- 



4. Harmonically Trapped Dipolar Bose Gases 

In this section we discuss the case of a harmonically trapped dipolar Bose gas, i.e., 
particles under the influence of the potential 

t/tr(x) = y {ujy + ujy + uy) , (so) 

where cuj denotes the trapping frequency in the z-direction. Due to the spatial 
dependence of the trapping potential (130|1 . the system is no longer translationally 
invariant and momentum is not a good quantum number anymore. Nonetheless, by 
means of the semiclassical and local density approximations (LDA) we will be able to 
derive analytical expressions for the physical quantities of interest such as, for instance, 
the condensate depletion and the equation of state. On top of that, we will investigate 
the influence of the quantum fluctuations upon the equilibrium configuration and the 
low-lying excitations of the system. 



4-1. Semiclassical and Local Density Approximations 

Let us start by implementing the semiclassical approximation to the non-local BdG 
theory. This can be done through the substitutions [HI HE] 

e^^£(x,k), ^ W (x, k) e*^ '', ^ V (x, k) e^^ '', (31) 

where the functions U (x, k) and V (x, k) are slowly varying functions of the position x. 
Within this semiclassical approximation, the BdG equations f lTSj) become 



e(x,k)- 



£(x,k)4 



2M 
~2M 



W(x,k)= v/M^y"d'xVint(x-xO v/^^ [V(x;k)+W(x;k)] e^<^'-^\ 
V(x,k) = v/M^y"dVVi,t(x-xO [V(x',k)+W(x',k)] e^''-^ 



(32) 
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The next step in order to solve the BdG equations fl32l) is to use the LDA for deriving a 
local term for the non-local dipolar interaction between the condensate and the excited 
particles. Denoting either Bogoliubov amplitude lA (x, k) or V (x, k) by g(x, k), the 
non-local term can be written in the semiclassical approximation according to 

/nl(x, k) = v^M^ j ^^x'VU^ - ^')^/M^q (x', k) e*k-(x-x'), (33) 

Under the LDA, this term reduces to 

JNL(x,k) ^e(x,k)g(x,k), (34) 
together with the abbreviation 

e (x, k) = gnoix) [1 + edd (3 cos^ ^ - l)] . (35) 

It is important to point out that the semiclassical procedure applied here can be 
justified within a systematic gradient expansion in the Wigner representation, where the 
local density approximation is shown to be the leading contribution [18]. This has the 
important consequence that systematic quantum corrections to the leading semiclassical 
term can be implemented, as in the case of the Thomas-Fermi model for heavy atoms 
[19]. Moreover, it also allows for estimating the range of validity of the LDA as explained 
in the appendix. 

By means of the LDA, the BdG equations ( fT5l) become simple algebraic ones, as 
in the homogeneous case. Therefore, the Bogoliubov spectrum can be obtained in the 
usual way and reads 

(x, k) = el^^ (x, k) - e (x, k) , (36) 

together with the definition of the LDA energy 

^LDA k) = ^ + e (x, k) . (37) 
Moreover, the semiclassical Bogoliubov amplitudes are given by 

" gLDA (X, k) 1 

r^i^ ' '\ • ^''^ 

We can now explore the effects of quantum fluctuations on interesting physical 
quantities such as the Bogoliubov depletion, the corrections to the ground-state energy, 
and the chemical potential. 



V^(x,k) = l 



4-2. Condensate Depletion, Ground-state Energy, and Equation of State 
Under the LDA, the depletion density reads 

, , , 8 / ?l(x)3a^ ^ / ^ , ^ 

An(x) = -^^-^Q3(edd). (39) 

Notice that the DDI enters this equation both directly, through the function Q3(edd), 
and indirectly, as it also determines the gas density n(x). Let us assume that the gas 
takes the shape of an inverted parabola with Thomas- Fermi radii Rx, Ry, and Rz, which 
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represents the mean- field solution in the Thomas- Fermi regime [HI [T2j, i.e., that the 
gas density is given by 



n(x) = n(0) 






(40) 



wherever this expression is positive and vanishes otherwise. In this case, the total 
depletion reads 

^ = '-^Vm^sQsie,,). (41) 

From Eq. fHT]) one identifies the gas parameter at the trap center, i.e. n(0)a^, as being 
decisive for the observation of the depletion in Bose gases. As for the dipolar contribution 
to the depletion, one should notice that increasing e^d up to the limit of stability of the 
ground state might increase the condensate depletion by only about 30%, see Fig [H 
Thus, any experimental observation seems quite difficult. 

Better prospects for an experimental observation of beyond mean-field effects are 
provided by the dipolar quantum correction to the ground-state energy. Indeed, the 
dipolar dependence of the energy density correction 



Ai?(x) = ^^7n(x)^/^^Q5(edd) (42) 

is controlled directly by the function Q5(edd), see FiglH For the sake of completeness, 
we shall also present the total correction for a parabolic condensate 

AE = b^gniO)^/MMQ5{e6d). (43) 

o 

In the following we shall see that the energy correction (1431) can be used for studying 
both the static and the dynamic properties of the system beyond the mean-field 
approximation. Equivalently, one can also use the quantum corrected equation of state 
of a dipolar Bose gas for this purpose. It can be obtained by differentiating the energy 
density with respect to the number density and reads 



/i(x) = UU^) + gn{^) + <|.dd(x) + ^n(x) W^^Q5(edd). (44) 

O V TT 

This equation obviously reduces to (12^ in the case of a homogeneous system. It also 
shows that there is no ambiguity in the chemical potential of a trapped system due to 
lack of translational invariance. Nonetheless, as in the case of a homogeneous system, 
only the mean-field contribution is anisotropic owing to the dipolar potential $dd(x), 
whereas the quantum correction, given by the last term in (jH]), is isotropic. The 
quantum correction remaining isotropic for a trapped system is an artifact of the LDA, 
as is explained in detail in the appendix. 
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4-.3. Variational Approach to Dipolar Superfluid Hydrodynamics 

Superfluity is characterized by the existence of an order parameter \1/ = ^/ne^^ whose 
modulus accounts for the superfluid density n and whose phase x accounts for the 
superfluid velocity. Therefore, it is possible to write down an action in the form [50] 



A[n,x] = - J dtd'^xnS^M x + ^^X^ +e[n]| 



(45) 



if one identifies the velocity with the gradient of the phase x according to v = Vx- In 
our case, the energy density e [n] is composed of a mean-field energy density 

cmf = f/tr(x) + |n(x, + dVKid(x - x')n(x', t) (46) 

and a quantum correction 

64 , , niy.,t)al _ . , 
eq = —gn{^,t)\ ^ ' ^ Qs ^dd • 47 
15 V vr 

As a matter of fact, extremizing the action (145|) with respect to the phase x(x, t) and the 
density n(x, t) leads to the continuity and the Euler equations (l27j) and ( 1281) . respectively. 
Therefore, this action contains all the elements that one needs in order to investigate the 
hydrodynamic properties of the system. However, there is a simpler and more efficient 
way of performing these studies than solving the aforementioned equations. Indeed, 
by choosing a special ansatz for the superfluid phase and density, one can address the 
physical properties of interest. This technique, slightly modified to include the Fock 
exchange term, has been applied before to study the hydrodynamics of dipolar Fermi 
gases both in cylinder-symmetric [51] and in triaxial traps [52] . 

We proceed with the extremization of the action by adopting a harmonic ansatz 
for the velocity potential 

X(x, t) = ^a,{t)x^ + ^ay{t)y^ + ]^a,{t)z\ (48) 

where the parameter controls the expansion velocity of the cloud in the z-th direction. 
Moreover, we use an inverted parabola as an ansatz for the particle density, which is 
given by 



n(x,t) = rao(t) 



' ^ x"^ y^ 



(49) 



Rlit) R^M\ ' 

wherever the right-hand side is positive and vanishes otherwise. Due to normalization, 
the Thomas- Fermi radii are related to the quantity fio(t) through 

15A^ . ^ 

Mt) = 50 

871 R{t) 

3 

with the geometrical mean R = R^RyRz- In order to render the notation more concise, 
the arguments of the functions are sometimes omitted, as long as no confusion can arise. 

By inserting the ansatz (jlHD and fH^ into the action (I45p . one obtains the action as 
a function of the variational parameters which allows to derive the corresponding Euler- 
Lagrange equations of motion. First, we obtain the equations for the phase parameters 
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to be ai{t) = Ri(t)/Ri{t). Then, with their help, we derive the equations of motion for 
the Thomas- Fermi radii of a dipolar Bose gas beyond the mean-field approximation. In 
the general case of a triaxial trap the equation for the Thomas-Fermi radius in the i-th 
direction reads 

Ri = —ufRi-\ I di{Rx, Ry, Rz, ^dd) + /J \ ■ (51) 

AuMRiR I R^^ J 

Here, we have introduced the abbreviation 

di{R,,Ry,Rz,edd) = l-edd[l-R^^R^]f (^^,^^ . (52) 
It includes the anisotropy of the DDI as expressed through the function 

/(..y)^l + 3. /,;^-''>-^f . (53) 

(1 — y^jVi- — a; 

where F {ip, k) and E {if, k) are the elliptic integrals of the first and second kind, 
respectively, with the arguments fc^ = (1 — ?/^)/(l — x^) and ip = arcsin \/l — x"^ |13]. 
It is worth remarking that the representation ( 153|) is valid for < x < y < 1. For 
other regions of the Cartesian plane, it has to be analytically continued |53]. For more 
information on the anisotropy function fl53|) . see Refs. [21 [511 

The quantum fluctuations are accounted for by the last term on Eq. ( ISTl) and their 
influence is characterized by the function 

/3(edd) = lQ,{edd)a'J'N'/', (54) 
where the numerical constant 7 reads 



/33 . 53 . 72 

7 = y ^ 4.49. (55) 

In their absence, the mean-field triaxial equations of motion of Ref. [2], which were first 
derived in cylinder symmetric form [TT], are recovered from fl3T]) . 

The beyond mean-field equations of motion (I5T!) represent the main result of the 
present paper. They allow us to investigate the effects of quantum fluctuations in a 
dipolar condensate in a triaxial harmonic trap. Indeed, solving these equations exactly 
is both difficult and unnecessary, due to the fact that the quantum corrections only have 
the particular form presented here, if they are small. For this reason, we will treat in the 
following all the /3 terms as small and calculate the physical quantities perturbatively 
only up to first order in (3. 



5. Cylinder-Symmetric Trapping Potential 

In practice, most experiments are carried out in traps which can be taken as cylinder 
symmetric to a very good approximation. For this reason, it is important to study this 
case carefully. To this end, we notice that the symmetry of the problem yields Ry = Rx 
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and we have to take into account the properties of the anisotropy function f l53p in the 
particular case /(x, x) = fs{x), which is given by |T2l [521 155] 

, , , 1 + 2x2 tanh"^ Vl - 

^^^^) = ^T^i^T^ ■ 

We remark that (156!) is vahd for < x < 1 and must be analytically continued for x > 1. 
The corresponding equations of motion (ISTj) . in this case, reduce to 



R.- -^.r^ + ^^^mWr 

Rz = — ^tR-z + 



2 15^iV 
" " ^ A-kMRIRI 



where we have introduced the auxiliary functions A and B. They depend only on the 
aspect ratio Rx/Rz and read 

5. i. Static Properties 

Let us first consider the effects of the beyond mean-field corrections on the stability of 
the system. It has been shown some time ago by means of a thorough mean-field analysis 
that a stable ground-state only exists for trapped dipolar condensates if the value of the 
relative interaction strength lies within the range < ead < 1 P=2j. For values of ead 
larger than 1, the ground state is, at best, metastable. Quantum fluctuations cannot 
alter this as their calculation within the Bogoliubov theory amounts to performing an 
expansion up to second order in the fluctuations of the field operators around their 
mean-field values. Such an expansion can only be carried out if the corresponding 
ground-state is stable. Therefore, the effects of quantum fluctuations on the properties 
of a dipolar condensate are only physically meaningful as long as < edd < 1- And 
this is clearly pointed out by the fact that both functions Q3(edd) and Q5(edd) become 
imaginary for edd > 1- 

As we discussed above, the most appropriate manner to study the effects of beyond 
mean-field corrections in equations (l57j) is to perform an expansion around the mean- 
field solution in powers of (5. To that end, we adopt the ansatz 

R, = Rl + SR„ Rz = Rl + SRz, (59) 

with being the mean- field Thomas- Fermi radius and 6Ri being a correction of the 
order j3, which is obtained by solving the static versions of (l57Il to first order in /3. Since 
all the quantities involved are functions of the Thomas- Fermi radii, the static properties 
of the system can be investigated by evaluating the corresponding correction. Consider, 
for example, the beyond mean-field aspect ratio 

K = ^ = K° (1 + 5k) • (60) 

R-, 
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Figure 3. Relative correction to 
Cdd for A — 0.50 (green), A = 0.75 
for Cdd = 0.89 (green), edd = 0.95 
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the aspect ratio Sk in units of Sk a) as a function of 
(red), and A = 1.00 (blue) and b) as a function of A 
(blue), and edd = 0.97 (red). 



Using ( l57l) with the Ihs set to zero, one first obtains for the mean-field aspect ratio the 
following transcedental equation [12] 



Then, by proceeding in the same way with (1571) up to first order in /3, the beyond 
mean-field correction of the aspect ratio is found 



5k 



2+eM[2{2-R,dRjB-{2 + R,dRjA]-2e^^^[A{l-R,dRjB + B{l + R,dRjA] ^^^o' 

(62) 

where the rhs is a function of the aspect ratio evaluated at its mean-field value. 
Moreover, we have introduced the quantity 

105v^ 



6K = ^^^aM0), (63) 

which sets the scale for the correction 6k of the aspect ratio. 

From fl62|) one recognizes that in the case of a non-dipolar Bose gas, for which k° = A 
holds according to ( 16T1) . the aspect ratio is not altered by the quantum corrections: 
though both Thomas-Fermi radii R^ and Rz are affected by the quantum corrections, 
owing to the isotropy of the contact interaction, one has SRx/R'^ = SRz/R^^y so that 
5k vanishes. In Fig. 13^), we plot the same correction as a function on the relative 
interaction strength edd at a fixed trap aspect ratio A. The green curve is for A = 0.50, 
the red one for A = 0.75, and the blue one for A = 1.00. For a vanishing dipolar 
interaction, the condensate aspect ratio is not affected by the quantum fluctuations, as 
we have discussed above. For increasing relative interaction strength edd, however, a 
non- vanishing correction shows up. When approaching the critical value edd = 1, above 
which the correction to the ground-state energy due to quantum fluctuations becomes 
imaginary, the correction to the aspect ratio remains finite though very large. This is a 
signal of the system becoming unstable. 
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Due to the presence of the dipole-dipole interaction, the role played by the trap 
anisotropy becomes an important feature of the aspect ratio correction. To exemplify 
this, we show in Fig. [3b) the correction to the gas aspect ratio 5k/ 5k as a function of the 
trap aspect ratio A for different values of the relative interaction strength eaa- Notice 
that the effect is larger for prolate (cigar like) traps and becomes smaller and smaller 
for oblate (pancake like) traps. 

In order to estimate the importance of the quantum correction to the aspect ratio, 
let us adopt the experimental values of the average trap frequencies and number of 
condensed particles from the ^^Cr-experiment reported in Ref. [3]. In that case, the gas 
parameter at the center of the trap is such that the unit of the variation of the aspect 
ratio is 5k ~ 0.05. This renders the observation quite difficult, as the aspect ratio 
variation would only become appreciable at large values of edd- For stronger magnetic 
systems, the situation is different. The s-wave scattering length of dysprosium, for 
example, is presently under investigation and there is evidence that it could be smaller 
than add = 133 cq [IE]. In this case, one would have edd > 1 and the present theory 
could not be applied, as this configuration would be, at best, metastable. Suppose 
that, by means of a Feshbach resonance, the scattering length of dysprosium could be 
set to Og = 150 ao- Then, one would have edd ~ 0.89 and, for the same number of 
particles and trap frequencies of the recently achieved dysprosium BEC [16] one would 
have 5k ^ 0.11, leading to much better prospects for observing these beyond mean-field 
corrections. 

5.2. Hydrodynamic Excitations 

In this section, we address the question of how the presence of dipolar interactions 
modifies the impact of quantum fluctuations on the low-lying excitations of a Bose gas. 

We proceed to calculate the shift in the excitation frequencies due to quantum 
fluctuations by separating each of the three Thomas-Fermi radii as a function of time 
in two contributions 

R^{t) = Ri{'d)+r]iSin{m + ^), (64) 

where -Ri(O) is the equilibrium value of the radius, rji represents a small oscillation 
amplitude of oscillation and VL the oscillation frequency. In addition, ip denotes a phase 
which is determined by the initial conditions. Notice that, instead of using from the 
cylinder symmetric equations ( 1571) . we actually go back to the triaxial equations (l5Ti) 
and, later on, evaluate the cylinder-symmetric limit. This procedure is necessary in order 
to study the radial quadrupole mode in addition to the monopole and the quadrupole 
modes. Then, one arrives at the eigenvalue problem 

^Oijirij = n^rii. (65) 

i 

In general, the matrix elements satisfy Oij = Oji. Moreover, since we are only interested 
in the cylinder-symmetric limit, we also have O^x = Oyy and O^z = Oyz- Thus, we are 
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left only with the following four independent matrix elements 

IbgN P RxdnJ. 



O 



xy 



hm col 

y^x 



lim col 

y^x 



lim 



lim 

y^x 



3 + 



IhgN (3 



AnMcol 



2RIR 



■9/2 



Ry 

Rx 
XR. 



+ 



15gN 



/3 



4 + 

Rzdn, dz 
d, + (3/R^'\ 
RxdRyd. 



^^Mul 2R^RyR"' 



dx 



IbgN 



/3 



According to prescription f l59|) . the matrix elements Oij in fl65|) are also corrected by 
terms of order (5. Therefore, we write them as 

0ij = 0l + 50i^, (67) 

with 50ij oc /3. To obtain the corrected oscillation frequencies we proceed as before 
and treat the terms of the order /3 as a perturbation. Expanding the corresponding 
frequencies up to first order in that term leads, at first, to corrected oscillation 
frequencies in the form 

VL = n^{l + 5VL) . (68) 

Requiring (!65|) to have non-trivial solutions leads to three values of the eigenvalue VL^ 
which correspond to the radial quadrupole, the quadrupole, and the monopole oscillation 
frequencies. In the following, we study the oscillation modes and discuss the perspectives 
for observing their corrections with respect to the mean-field values. 

5.2.1. Radial Quadrupole Mode Let us consider, at first, the eigenvalue of f l65|) which 
corresponds to the radial quadrupole mode. For this mode, one has rj^ = —% and 
r]z = 0. The oscillation frequency of this mode is given according to 

(69) 

one obtains the mean-field radial 



-P/r'^\ 

RxdR,dx 

^3/2 



(66) 



rq 



^/0 XX '^xy 

By evaluating these matrix elements from 
quadrupole frequency which can be written as 



rq 



CO, 



2 + edd/irq(K' 



with the abbreviation 



/irq(K) 



3 ,[2(1 



-K 

4 



(4 



(1 - t^^Y [1 - eddA(/s:)] 



(70) 



(71) 



Previous studies of the radial quadrupole mode of dipolar Bose gases have been 
carried out perturbatively [56] and also at the mean- field level [57]. Notice that in the 
case of a pure contact interaction, i.e., for edd = 0, the mean- field radial quadrupole 
frequency does not depend on the geometry of the system at all. Correspondingly, a 
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Figure 4. Quantum correction to the frequencies of the low-lying excitations in 
units of (5fi as a function of the trap aspect ratio A. The colors black, red, and blue 
correspond, in turn, to the nionopole, quadrupole, and radial quadrupole oscillations. 
Moreover, continuous curves are for the parameter values = lOOao and e^d = 0.69, 
whereas dashed curves represent edd = [55]. In the latter case, the curves do not 
depend on a^. 



quantum correction can only have its origin at the presence of the DDI. Indeed, the 
anisotropy of the dipolar interaction is responsible for modifying the radial quadrupole 
frequency according to the expression 



2f]i OK ri _^,A(^.\^ m^' 



(72) 

I [1 - e^^Ain)] ■ " 

which immediately vanishes for non-dipolar Bose gases. This can also be clearly seen 
from Fig. HJ where the correction SQ^q is shown in units of 

S~^ = ^^/^^W) (73) 
as a function of the trap aspect ratio A. Here, we consider, for instance, erbium, which 
has a magnetic dipole moment of m = 7yUB and assume it to have an s-wave scattering 
length about the same as in ^^Cr, i.e., = lOOao, yielding edd ~ 0.69. Moreover, we 
have adopted realistic values for the particle number and the average trap frequency 
from Ref. |3], for which one obtains SQ ~ 1%. In the absence of the DDI (dashed 
curve), the quantum correction to the frequency vanishes, whereas it is non-zero in its 
presence (solid curve) and might amount up to 1 or 2 percent. This represents a clear 
signal for detecting many-body effects stemming from the DDI in cold atomic systems. 

5.2.2. Monopole and Quadrupole Modes Let us now turn our attention to the other 
two modes. They correspond to oscillations in which the radial and axial coordinates 
vibrate either in phase, as in the case of the monopole mode {t]^ = Vy ^ Vz), or out 
of phase, as in the case of the quadrupole mode {rj^ = % ^ ~Vz)- Therefore, these 
modes are denoted with a plus and a minus index, respectively. In accordance with the 
previous reasoning, we write the oscillation frequencies in the form 

n± = (1 + 6n±) , (74) 
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where denote the exact mean-field monopole (+) and quadrupole (— ) frequencies 
and 5VL± denotes the relative quantum correction of order /3. 

The mean-field frequencies have been investigated by O'Dell et al. [H]. Adapting 
our triaxial notation to their cylinder-symmetric one, the mean-field frequencies can be 
written as 



(75) 



with the mean-field matrix elements 



1 + edd 



0° +0° 

^ XX ' xy 



col + 3ul- 



2, .2 



k2( 1+4^2 
2{1-k2)2 



Cdd 



1 - 

5~2k2 



_ 3(4+>c2)/4k) 
2(1-k2)2 



eddA{K) 



10^ 



2h' 



20 



2ujIk- 



1 + edd 



1+2^2 I 15k2/3(k) 
■ 1_k2 -1- ..9x2 



2(1-k2)-= 



1 - edd^(/«) 

We now take advantage of the fact that the mean-field matrix elements are functions of 
the aspect ratio k = R^/Rz alone and not of the radii individually. This allows us to 
calculate the contribution to the corrected eigenvalue problem due to the change in the 
aspect ratio. In addition, there is a further contribution coming from the fact that the 
equations of motion have themselves been corrected. Together, both contributions are 
given by 



^hxx 
Shzz 



edd {I - RzBr^) A] 



Ok 



6k, 



6k + 



/3k 



(1 - edd^)' 
2[l-edd {5A + 8B-RzdRM 



2RlRf'^ 
(5ujI [1 



edd (1 



;i-eddA)' 
-2RzdR^)Az 



2RW^'^ 



[I - ed^AY 



Finally, for the relative correction to the frequencies, we obtain 



6VL 



± 



6Kx+SKz ± 



2 {K,jKx+KJKz) + {Kx-hl) {6Kx-6K.) 



(77) 



(7J 



^^xz^^x~'^i^xx ^^zz) 

In order to appreciate the effect of the quantum corrections on realistic experimental 
systems, we plot in Fig. H] the corrections of the monopole and quadrupole frequencies 
as functions of the trap aspect ratio A in units of 6VL. As we remarked in the discussion 
about the radial quadrupole mode, typical experiments have 6VL ^ 1%. Thus, for 
example, the quantum correction for the monopole oscillation frequency (shown in black) 
of a moderate cigar shaped trapped gas could amount to as much as 5%, so that one can 
realistically expect this effect to be measurable. In Fig. HI the dashed lines correspond 
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Figure 5. Thomas-Fcrmi radii (blue) and Rz (red) in arbitrary units as functions 
of time in units of o;^^ during monopole oscillations. Wo adopt the values edd = 0.89 
and as ~ lOOao- In part a) we show the oscillations from u^t ~ until uj^t = 5, while 
in part b) the plots go from uj^t = 40 to uj-^t = 45. 

to non-dipolar Bose gases, i.e., to edd = 0. It is interesting to observe that the presence 
of the DDI changes these curves qualitatively. For the monopole and quadrupole curves, 
the DDI leads to a crossing of the corrections at some value of A, which depends on 
the relative interaction strength edd- The fact that, for given values of edd and A the 
quantum correction of the monopole frequency becomes smaller than the correction of 
the quadrupole frequency is absent for non-dipolar Bose gases and represents, therefore, 
a clear signature of the DDI. We remark that this feature is present even for weakly 
dipolar systems such as chromium. 

In view of the recent important experiment, in which Bose- Einstein condensation 
of dysprosium was achieved [I6], we plot the Thomas- Fermi radii as a function of time 
for monopole oscillations in Fig. [5] and for quadrupole oscillations in Fig. |6l Thereby, we 
have adopted the values of the number of particles and trap frequencies from Ref . [16] , 
which give a trap aspect ratio A = 3.8. Moreover, by choosing the s-wave scattering 
length of = 150ao, we have a relative interaction strength of edd = 0.89. As a matter 
of fact, the amplitudes if the oscillation for and Rz are not independent from each 
other. Instead, their ratio is an intricate function of the system parameters. However, 
this is irrelevant for the effect that we aim for and we take each of the amplitudes to be 
10% of the corresponding radius, i.e., r/j ^ 10%i?°. 

To analyze the monopole oscillations, consider the Thomas-Fermi radii R^ and R^ 
as functions of time from u^t = to co^t = 5, shown in Fig. |5^), and from u^t = 40 
to u^t = 45, which appear in Fig. [5]d). The only difference, which is initially visible, is 
the fact that Rx becomes larger due to quantum fluctuations. For R^, for example, the 
correction is too small to be seen here. Moreover, the valleys and hills do seem to match 
very well. In the course of time, however, the mean-field and the quantum corrected 
oscillations do depart from each other. In the case of the quadrupole oscillation, which 
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Figure 6. Thomas- Fermi radii Rx and Rz in arbitrary units as functions of time 
in units of uj^^ during quadrupole oscillations. We adopt the values edd = 0.89 and 
as ~ lOOflQ. 

is depicted in Fig. [6l the curves can be distinguished from each other at much smaller 
times even for R^, which possesses very small amplitude corrections. The reason for 
that lies in the values of the trap aspect ratio A = 3.8. As one can see in Fig. HJ for this 
particular value of A, the correction of the quadrupole frequency is much larger than the 
one of the monopole frequency. Thus, as we have just demonstrated, in order to achieve 
the detection of many-body effects on the low-lying excitations of dipolar Bose gases, it 
might be more adequate to consider monopole or quadrupole oscillations depending on 
the trap aspect ratio A. 

6. Conclusions 

We have theoretically investigated beyond mean-field properties of both homogeneous 
and harmonically trapped dipolar Bose gases, focusing on the low-lying excitations. 
After having studied the Bogoliubov-de Gennes theory, we have characterized the 
influence of the DDI on the condensate depletion, on the equation of state, and on 
the sound velocity of homogeneous Bose gases. With the help of the local density 
approximation, these results could be generalized to the case of a harmonically trapped 
gas. Then, within the framework of superfluid hydrodynamics, we have variationally 
derived equations of motion for the Thomas-Fermi radii and used them to investigate 
the case of a cylinder symmetric trap. While difficulties in performing precision 
measurements of the particle density represent a hurdle for identifying dipolar beyond- 
mean field effects in static properties, the oscillation frequencies offer much better 
perspectives. The radial quadrupole mode, for example, acquires a finite quantum 
correction which clearly has its origins in the DDI. The frequencies of the other two 
modes are also modified both quantitatively and qualitatively by the inclusion of the 
DDI in the beyond mean-field regime. As a result, the low-lying oscillations of Bose 
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gases offer an important possibility for observing many-body dipolar physics. 

The beyond mean-field theory in the form presented here could be adapted to study 
a variety of important problems. An obvious suggestion is a systematic investigation 
including, e.g., the scissors mode beyond mean- field approximation [57]. Moreover, 
presently available studies of the dipolar dirty boson problem, for example, concentrate 
on homogeneous systems [59] and on the equilibrium in harmonic traps [60] , but nothing 
is yet known about the dynamical aspects of the trapped system. Also, the non-linear 
dynamics induced by means of modulation of the s-wave scattering length [61] would for 
sure have a non-trivial interplay with dipolar interactions which is still to be investigated. 
In addition, it would be interesting to clarify how quantum fluctuations can alter the 
character of dipolar systems with the dipoles not aligned along the main axis [62] or set 
to rotate [63l . 
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Appendix: Validity of the Local Density Approximation 



In this appendix we will clarify under which circumstances the LDA can be used for 
long-range interactions such as the DDI. To this end, we will use the next order of 
the corresponding gradient expansion in order to estimate the error brought about by 
neglecting terms of higher order than the LDA term. 
Consider the non-local term in Eq. ( ITSll 



/nl(x, u) 



'rio^xj 



rao^x'jg,, (X 



(A.l) 



where (x) stands for either Bogoliubov amplitude lAy (x) or Vy (x), and recall that v is 
a discrete quantum number while no(x) denotes the condensate density. In this case, the 
semiclassical approximation is obtained under the substitution ( 13T|1 with g(x, k) being 
a continuous and slowly varying function of x and k. Then, /nl(x, k) becomes 

Jnl(x, k) = j d=^x'Mnt(x - x')F(x', k)e^^ (^-^'), (A.2) 

^no(x)g (x, k). Performing the variable transformation x' — !■ x' + x. 



with F(x, k) 
one has 



/NL(x,k) = v/^o(x)<^F(x,k) 



+ 



Vlnt(k) 



(A.3) 



where the gradients act only inside the curly brackets and in the direction of the arrows. 
Moreover, the dots replace higher order terms which are neglected consistently with the 
Thomas-Fermi approximation. 



Beyond Mean-Field Low-Lying Excitations of Dipolar Bose Gases 



25 



For a contact interaction, we have VkMnt(k) = 0, so that the LDA is immediately 
justified. In the case of the DDI, however, this is not the case. First, notice that the 
leading order of the gradient expansion is isotropic, but the next-to-leading term is not, 
as the spatial variation of the density is coupled to the k space. This explains why the 
anisotropy of l^nt(k) does not directly affect the quantum corrections at the LDA level. 
Notice that, in the usual experimental case of large particle numbers in which we are 
interested, the next-to-leading term can be neglected. Indeed, its ratio to the LDA-term 
can be estimated by the substitution Vx l/i?TF, Vk V-^c, with i?TF the mean 
Thomas- Fermi radius and Kc the momentum scale defined by the speed of sound. These 
are more conveniently expressed in terms of the condensate density at the trap center: 



R 



TF 



^7^0(0) 

Then, one finds that the LDA is valid as long as 

clda 1 



M', = v'MV^„t(k)no(0). (A.4) 



< 1 (A.5) 



[N^aln,{Q)Y/^ ^1 + e^d (Scos^^- 

with the numerical constant clda = (3^5^7r)^^/^. 

As one would expect, for a large enough particle number iV, the LDA is valid as 
long as the system is stable, i.e., for relative interaction strengths satisfying e^d < 1- 
Indeed, the factor ci,da[N'^o?s^o{^)\~^^^ can be expressed as A/2(aho/15A^as)^/^ with the 
average oscillator length Cho- For typical experimental situations with non-dipolar gases 
the Thomas- Fermi parameter fulfills Nas/a^o S> 1, this being the condition of validity of 
the hydrodynamic description. Moreover, even in dipolar experiments with chromium 
(see Ref. [3]), in which the s-wave scattering length is decreased on purpose in order 
to enhance the relative dipolar interaction strength, one might have Nas/a\^Q ~ 60 for 
the smallest particle numbers and scattering lengths. Thus, we conclude that the LDA 
is an applicable approximation in many situations of interest for the physics of dipolar 
Bose gases. 
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